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Abstract 

We consider a variant of hybrid inflation where the waterfall phase transition happens during inflation. 
By adjusting the parameters associated with the mass of the waterfall field we arrange that the phase 
transition is not sharp so inflation can proceed for about 50-60 e-folds after the waterfall phase transition. 
We show that one can work in the limit where the quantum back-reactions are subdominant compared to 
the classical back-reactions. It is shown that significant amount of large scale curvature perturbations are 
induced from the entropy perturbations. The curvature perturbations spectral index is either blue or red 
depending on whether the mode of interest leaves the horizon before the phase transition or after the phase 
transition. This can have interesting observational consequences on CMB. The non-Gaussianity parameter 
fNL is calculated to be < 1 but much bigger than the slow-roll parameters. 
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INTRODUCTION 



Recent observations [jj strongly support inflation as a correct theory of early universe and 
structure formation [2]. Thanks to the precision data different inflationary models can be distin- 
guished based on their predictions for the curvature perturbation spectral index, the amplitude of 
the primordial gravitational waves and the level of non-Gaussianity. 

It proved very difficult to obtain an appreciable amount of non-Gaussianity in simple models of 
inflation. One has to have either multiple field inflationary scenarios or non-trivial sound speed, 
for a review see [4], Q], Q], Q], [gj], ij] and the references therein. Specifically, in models of 
multiple field inflation when the slow-roll conditions are violated temporarily on the field space, 
one may naively expect that an appreciable amount of non-Gaussianities can be produced. Careful 
examinations in the context of double inflation indicate that this may not be the case [lo| . 
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Therefore, it would be interesting to extend these analysis to similar models where non-trivial 



dynamics such as a sudden change in sound speed during inflation, fields annihilations [12. 



particle creations 



end of inflation 
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17j and phase transition 18|-2(| happening either during inflation or at the 



27]. 



281 ] , [29j where the waterfall phase transi- 



In this work we consider a variant of hybrid inflation 
tion happens during early stage of inflation. By tuning the effective mass of the waterfall field, we 
arrange that the waterfall phase transition is mild enough such that inflation continues for a long 
period, say 55 e-foldings 



3d . l3ll ] . To bring the effects of phase transition into cosmic microwave 
background (CMB) observational window, we shall assume that the phase transition happens 
around first few e-foldings, say first five e-folds. We follow the dynamics of both fields so our 
treatment is a two-field inflationary mechanism throughout. We would like to see, first, whether 
the entropy perturbations can induce significant amounts of large scale curvature perturbations, 
and second, whether a significant amount of non-Gaussianity can be produced. These questions 



32, 



33] attracted new interest in the literature for the model of standard hybrid inflation where the 



waterfall phase transition happens very efficiently at the end of inflation. As demonstrated in 

n 
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351 ]. [361 ] and [37J] the quantum back-reactions from very small scales inhomogeneities produced 
during the waterfall phase transition uplift the tachyonic instability and shuts off inflation very 
efficiently. As a consequence, the large scale curvature perturbations are exponentially suppressed. 
By the same reasoning, it seems natural to ask whether this conclusion can be averted if one relaxes 
the model parameters such that the waterfall phase transition, happening during the early stage 



of inflation, is not very sharp. This is one of our main goal in this work which we shall elaborate 
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in details in the subsequent sections. 

As a remark, somewhat related to our work here, there have been many works in the literature 
concerning curvature perturbations and obtaining non-Gaussianities in models where there are 
local feature during inflation. In 



38, 



condition in a single field model. In 



39 



40, 



this was translated into a sudden violation of slow roll 



4l[ the starting model is a multiple field scenario, where 



the power spectrum features are induced from slow-roll violation or phase transition due to other 
non-inflationary fields. However, the analysis of curvature perturbations and non-Gaussianities in 
these models are as in single field inflationary models. 

The rest of the paper is organized as follows. In section [II] we present our background and the 
classical fields equations. In section IIIII we study the quantum excitations of the waterfall field in 
details and calculate the power spectrum of the entropy perturbations. In section HVl we compare 
the quantum back-reactions and the classical back-reactions and specify the limit where the former 
can be safely ignored compared to the latter one. In section|V]we study the curvature perturbations 
induced from the entropy perturbations in details. In section IVII we use the complementary 8N 
formalism to calculate the level of non-Gaussianity in our model. Conclusion and discussions are 
given in section IVIII followed by an appendix describing limiting behaviors of the Bessel functions 
used extensively in our analysis. 



II. WATERFALL PHASE TRANSITION DURING INFLATION 

As explained above we consider a variant of hybrid inflation model where the waterfall phase 
transition takes place during inflation. The potential is 

v(<p, ^) = \{^-^Pj + \™ 2 <? + ^ W , (i) 

where (j> is the conventional inflaton field, tp is the waterfall field and A and g are two dimensionless 
couplings. The system has a global minimum given by <p = and ip = M/yX. Our assumption 
is that the first period of inflation takes place during <p c < <f) < where (f>i is the initial value 
of the inflaton field and 4> c = M/g is the critical value of 4> where the waterfall field becomes 
instantaneously massless. Unlike conventional hybrid inflation, we assume that the waterfall field 
is not heavy so it can also slowly roll down during the first stage of inflation. For (j) < (f> c the 
waterfall becomes tachyonic triggering an instability in the system. Again, unlike standard hybrid 
inflation, we assume that this phase transition is mild enough such that it will take a long time 
for the waterfall field to settle down to its global minimum. This provides the second stage of 
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inflation. Our parameters should be such that the second stage of inflation is long enough say 55 
e-foldings or so. 

It is well-known that the potential JT]) with the discrete Z2 symmetry is plagued with the 
domain-wall formations at the end of inflation 42] which are dangerous cosmologically. To get rid 



of this problem one has to consider ip to be a complex scalar field so at the end of inflation cosmic 
strings are produced which can be safe cosmologically. However, the process of topological defect 
formation is beyond the scope of this work and we shall proceed with ip to be real as is customary 
in many models of hybrid inflation in the literature. 

A. The Background Fields Dynamics 

Here we study the classical evolutions of background fields <p and ip before, during and after the 
phase transition. The background space-time metric, as usual, is 

ds 2 = -dt 2 + a(i) 2 dx 2 , (2) 

where a(t) is the scale factor. 

It is more convenient to use the number of e-foldings, N, as the clock dN = H dt. The 
background fields equations are written as 

4>" + (3 - e)<f>' + (a + g 2 ^j = (3) 

iP" + (3 - e)iP' +(-f3 + g 2 ^ + X^\ ip = Q, (4) 

where the prime denotes the differentiation with respect to the number of e-foldings and e = 
—H/H 2 is the slow-roll parameter which is assumed to be very small by construction. Also the 
dimensionless parameters a and (5 are defined as 

m 2 m 2 3Xm 2 m 2 

a = w ' ao = ^ = W (5) 



and 



where m p = 1/G for G being the Newton constant and -f/o = \/2ir/3X M 2 /m p . Here the quantities 
ao and (3q respectively represent the values of the parameters a and f3 in the limit when we 
neglect the variation of H during inflation so H ~ Hq, a ~ ao and /3 ~ /3q. However, in order 
to get accurate enough solutions it is important to consider the running of H in our analytical 



treatments. We work in the limit where a <C 1 corresponding to a light eft field. Also since we 
are interested in a mild phase transition, we take /3 < 1. This is in contrast to standard hybrid 
inflation model with a sharp phase transition at the end of inflation where /3 — /3q S> 1. 

To simplify the notation, we take the critical point as the reference point and define n = N — N c . 
We use the convention that at the start of inflation for <f> = (pi, N = 0, at the time of phase transition 
N = N c and at the end of inflation N = N e . With this convention n < before the phase transition 
whereas n > afterwards. 

Let us for the moment assume that a and (3 are constants. At the end of analysis, we will include 
the running of a and /3 effectively in the analysis. Also we are in the limit where g 2 ip 2 /H 2 <C a, 
i.e. the back-reaction of the waterfall field on the inflaton field is small during inflation. With 
these assumptions one can easily solve Eq. Q to get 



with 



{n) ~ 4> c exp (-r n) (7) 



'•<>- I tj-\/t -''"-I -Y ( 8 ) 




Equivalently, one also has 

N c ^-\n(^]. (9) 

As explained before, one can not neglect the running of a and (3 which results in significant 
errors in the results above. Here we take into account the running of a and /3 . For this purpose, 
we consider the next leading term in H from the Friedmann equation 

H a ^ ( 4A + 2 m * ) ■ (10) 

It turns out that the correction to a are less crucial as compared to j3 and an overall averaging 
would suffice. We define the averaged H 2 via 

1 i-N 

H 2 = —J H 2 {n)dn. (11) 
Plugging the result from Eq. (|10p and performing the integral, one obtains 

Using this averaged value of H 2 into the definition of a results in 



in which N e is total number of e- foldings. This equation modifies ro to an effective value of r = a/3. 
Now plugging Eqs. (|7J) and (|12p in the definition of (3 (with ro — > r) one obtains the first order 
correction to (3 

P^P {l-Te- 2rn ) , (14) 

„ 4vr / C \ 2 

where F = — «o < 1. 

3 \m p ij 

Let us now turn to the dynamics of ip field. As it can be confirmed from our full numerical 
results we are in the limit where if) 2 /Hq <C /3/A so the equation for ij) simplifies to 

< + 3^ + /3(e- 2rri -l)^ = 0. (15) 

Using the modified form of (3 from Eq. (|14p . the equation of motion for ip modifies to 

V>" + 3v' - A) [i - (i + r)e" 2rn + r e - 4 ™] v = o . (16) 

The solutions of the above equation are given in terms of the Whittaker functions 
e (r-3/2)« M ^ (YNL e -2rr^ and e (r- 3 /2)n W ^ M^ e -2rn\ in which K = 1/ = 

— \fPo + 9/4 and ^ = v/2. In the solution above one can check that k> 1 and we can approxi- 

r 

mate the Whittaker functions with the Bessel functions and 



(17) 



$ ~ e- 3 ™/2 

where f) = /3 (1 + T). 

S _ 

For our range of parameters one can easily show that v > — — e rn 3> 1 and in this limit Y v (x) 

r 

is much larger than J u {x). After imposing the initial conditions the contributions from J v in Eq. 
(fTTj) becomes negligible and the classical trajectory can be approximated by 



e^/ 2 f . (18) 

Y v { &e* 



VI„rN c 

We have checked that this analytic formula for ip is in good agreement with the results obtained 
from the full numerical analysis. 

III. DYNAMIC OF QUANTUM FLUCTUATIONS 

In this section we study the dynamics of waterfall fields quantum fluctuations which play the 
role of entropy perturbations. The goal is to calculate the curvature perturbations induced from 
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these entropy perturbations which would be the subject of the studies in section [V] As we shall 
see the adiabatic curvature perturbation from the inflaton field is subdominant compared to the 
curvature perturbation induced from the entropy field. 

The equation governing the dynamics of waterfall quantum fluctuations, Sipk, in momentum 
space k is 

5^ + 35^+ (J^~P + g 2 ^\ 6^ = 0. (19) 
By substituting the equation of (j> field from Eq. [7] one has 

<ty£ + 3#k + + P ( e ~ 2rn ~ !)) = . (20) 

Since the effective mass of Sip^ is at the same order as H , that is j3 < 1, one can neglect the term 
containing f3 in Eq. (|20p for the sub-horizon perturbations and the solution of the <5Vk excitations 
inside the horizon is given in terms of the Hankel functions H^ 2 (ke~ n / k c ) and H^ 2 (ke~ n /k c ) 
where k c = He c is the critical mode which exits the horizon at the moment of phase transition. 
We require that deep inside the horizon the solutions start from the Bunch-Davis vacuum 

e -i/cr 

5ibh — > — t= as — kT — y oo , (21) 
k aV2k 



where t is the conformal time, dt = —adr. With this initial condition the incoming solution for 
the modes inside the horizon is obtained to be 



= -J-Zr *- N ° e~ 3n/2 H$> 2 ± ) . (22) 



As can be seen from this expression, the amplitude of the quantum fluctuations at the time of 
horizon crossing n* when k = e n "k c is given by 

|^|~-^=. (23) 
V Ak 6 

Note that here and below an asterisk represents the values of the corresponding quantities at the 
time of horizon crossing. 

After horizon crossing one can neglect the term containing k 2 in Eq. (|20p and the equation for 
V'k becomes identical to the background ip equation, Eq. (I15p . with the answer similar to Eq. [18] 



ciJ„ ^e~ rn +c 2 Y u 



(24) 



We need to fix the constants of integrations c\ and c 2 by imposing the matching conditions connect- 
ing the outgoing solution Sipt to the incoming solution Srpr . The matching condition is performed 



at n = n m when the term containing (3 in Eq. (|20p becomes comparable to the term containing k 2 
which results in the following equation for n m 



-2n„ 



f3\e 



-2rn„ 



II 



2(3 r\n rt 



(25) 



From the above equation one observes that n m can be either positive or negative for the physically 
relevant modes. This means that the time of matching can occur either before or after the waterfall 
phase transition. Furthermore, comparing Eq. (|25]) with the equation of n*, that is e n * = k/k c , 
one concludes that n m > re* and for a given mode the time of performing the matching condition 
is always after the time when that mode leaves the horizon. 

To further simplify the analysis of matching conditions we can neglect the running of n m com- 
pared to e™ m and replace Eq. (|25|) with the following simpler equation 



k 



V2A 



or. 



Imposing the conditions 6i/j k (n m ) = 6ipj~(n m ) and dtp. 



ci 



C2 



7T 7re 
k, 



-N c 



7T 7TB 



4r 

-N c 



-rn m rr(l) v ' 
-"3/2 l v 



(26) 



Sifjj^(n m ) we can fix c\ and C2 



""3/2 Iv 



'3/2 



k 
k r 



(27) 



Here primes denote derivatives with respect to the argument of the Bessel and the Hankel functions. 
Also the arguments of the Hankel and the Bessel functions, respectively, are the same as those in 
(|22p and (|24p with n = n m . We have checked that with these values of c\ and C2, our analytical 
solutions (|22p and (|24|) are in very good agreements with the results obtained from the full numerical 
analysis. 

In the Appendix we presented approximate formulae for c\ and C2. Considering the fact that 



we are in the limit where v > \J /3o(r + l)e rn /r 3> 1 and using the approximate expressions for 
ci and C2 one can check that the term containing J v in (I24p becomes subdominant and 



HI 



c 2 e 



-3n/2 v 



(28) 



Now it is time to compute the power spectrum of entropy perturbations, S. Following the 
prescription of 43f] we can perform a local rotation from <j> — ip field space into a — s field space 



where da = cos 6 d<p + sin 9dip represents the adiabatic field tangential to the classical trajectory 
while ds = cos 9dvp — smOdcp represents the entropy field orthogonal to the classical trajectory. 
Here 6 = tan -1 ^' /(f)') is the angle between (f> and ip in the field space [43|]. As can be verified 



from our numerical analysis we are in the limit where the trajectory in ip — (f> space is very flat, 
corresponding to 6,6' < 1. We will further elaborate on this point in section [V] This implies that 
<5sk ~ and = -f-^Sk ~ which, using ([28]) . results in 

P ^^Vl^! N2l 4£-»le-. (29, 



i I (2ir) ! 

The ^-dependence of the entropy power spectrum is only due to the pre-factor A; 3 and the constant 
of integration C2- Using (|26p one has din A; = dn m so the spectral index of entropy perturbation, 
n s , in terms of n m is 

n s -l = 3 + ^^. (30) 
dn m 

In order to compute the spectral index analytically we use the approximate expression for C2 given 
by Eq. (|A9|) an d the limiting behavior of Bessel function given by Eq. (jA7]) to obtain 

n s -l.f>(f . (31) 

We have checked that this gives qualitatively a good approximation for n s when compared to the 
full numerical analysis. As explained before, n m can be either positive or negative. For modes 
which leave the horizon before the phase transition n m < whereas for modes leaving the horizon 
after the phase transition n m > 0. This means that n s can change from a blue spectrum to a red 
spectrum, depending on whether the mode of interest leaves the horizon before the phase transition 
or after the phase transition. This conclusion has been verified numerically. 

IV. BACK-REACTIONS: QUANTUM OR CLASSICAL ? 

Before we proceed to calculate the curvature perturbation and its power spectrum, we have to 
determine whether or not the quantum back-reactions are small compared to the classical back- 
reactions. In standard hybrid inflation model with ^ > 1, corresponding to a very sharp phase 



transition, it was shown in 



34 



- 37j] that the quantum back-reactions from very small scales inhomo- 
geneities produced during the waterfall phase transition dominate exponentially over the classical 
back-reactions. The back-reactions of quantum fluctuations uplift the tachyonic instability during 
the waterfall phase transition and shuts off inflation very efficiently. In our case at hand with /3 < 1, 
corresponding to a mild phase transition, one may expect that depending on model parameters the 
quantum back-reactions are sub-leading and one can only use the classical back-reactions induced 
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from Aip 41 and g 2 ^ 2 (p 2 interactions to terminate inflation. As shown in [35|, the latter becomes 
important slightly sooner than the former. 

The expectation value of the quantum fluctuations sometime after the phase transition is 

W 2 ) = /g^. (32) 

As mentioned before 9 <C 1 during most of inflationary period and <5sk — 5^k- We already 
calculated for super-horizon modes, given by Eq. (|28p . We also note that for super- horizon 
modes 8ip^ evolve as the background tp field given by Eq. (|18|) such that 5sk(N) = i}(k)ip(N) 
where 



m^^YA^e^) e 3 ^ 2 (33) 
tpi \ r J 

Here we provide an approximation for Q(k) which helps us to evaluate the integral in Eq. (I32p . 
Using Eq. IA9I one has 

P2 7re~ Nc 

|C2l " 7(x) V^^ J ^ (x) ' (34) 



in which 



VA)(r + l) e _ r „ m>>1 ^ -n™ ~^2^<i ) (35) 

r k c 



and 



7(x) 



, 2i/ 2 - x 2 , 
rfz + z) + 3z + 3i 



~ 6* + C(/3 2 ) . (36) 



One can check that the main fe-dependence of C2 comes from the Bessel function and in our 
approximations 

m * ^=tAt * ■ ( 37 ) 

As a measure of the strength of the quantum back-reactions, we calculate the ratio — — and 

V ( n ) 

see under what conditions this ratio is small so one can safely neglect the quantum back-reactions. 
Using the above approximations one has 

m^Hi r k f d^^ = jii_ kf_ 

^ - $ 4 (2^)3 2fc3 4^2 k . ■ W 

Here k{ and kf, respectively, correspond to the largest and smallest modes which become tachyonic 
during inflation. For the smallest scale which becomes tachyonic during inflation we have kf ~ 

10 



y / 2(3rnf exp (Nf) H and for the largest mode we can set ki = H. Plugging these into Eq. (f38|) 
and noting that y / 2/3rnj < 1 results in 



(5^ 



(39) 



lp 2 l/jf 47T 2 ' 

This equation shows that the quantum back-reactions can be safely ignored if Hq < ipi, that is if 
one starts with large enough classical waterfall field values at the start of inflation. It would be 
more instructive to express the ratio Hi/ipf in terms of the number of e- foldings and the mass 
parameters. Using Eq. (|18p one obtains 

2 



9 



~ — e 

a 



-3JV> 



of the waterfall field on the inflaton field 



~ m 



(40) 



To get this relation, it was assumed that the end of inflation is determined by the back-reactions 

1 '~ Using the approximations for the 



35fj so g 2 i/j 2 



Bessel functions given in Eq. (|A7|h the ratio above is simplified to 



Hi 



9 2 

— — exp 

a 

9 2 

~ — exp 

a 



(2ris 
2/3 



3)N f 



„2rN c 



2vr 2 



N 



f 



1 

2r 



(41) 



Combined with Eq. (|39|) . the condition under which one can safely neglect the quantum back- 
reactions till the end of inflation is translated into 

-,2 



9 

— <C exp 

a 



2r 



exp [-30/3] , 



(42) 



where the final approximation is for typical values of r and Nf used in our numerical analysis. For 
our numerical example with /3 = 0.7, a ~ 0.04 and g 2 = 2 x 10~ 12 this condition can be met easily. 

Eq. (|42p indicates that the strength of the quantum back-reactions is exponentially sensitive to 
the parameter /3. For fixed values of the inflaton mass and coupling g, one has to start with small 
enough parameter (3 such that the quantum back-reactions can be safely ignored. This conclusion 
is consistent with our starting intuition that if the phase transition is mild enough one can neglect 
the quantum back-reactions. This is also consistent with the conclusion drawn in the model of 
standard hybrid inflation with /3 3> 1 that the quantum back-reactions dominate exponentially 
over the classical back-reactions 35]. 
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V. CURVATURE PERTURBATIONS 



We now have all the materials to calculate the final power spectrum of curvature perturbations. 
For this purpose we need to know the amplitude of adiabatic curvature perturbations at the time of 
phase transition as the initial condition and integrate the evolutions of curvature perturbation from 
the time of phase transition till end of inflation. The final amplitude of curvature perturbation, 
therefore, is 

n f =K + / K'dn, (43) 
Jo 

where TZo represents the adiabatic curvature perturbations in the absence of entropy perturbations. 
As demonstrated in [43J the evolution of curvature perturbations for the super-horizon modes, 



induced by the entropy perturbations, can be written as 

TZ' = ^-5s. (44) 
a 

As one can see from above equation both 6' and 5s can source the curvature perturbations. We 
also note that 6' represents the acceleration of ip, specially during the phase transition. As can be 
seen from our full numerical analysis, the classical background is such that during inflation and 
phase transition, 6,6' -C 1- Inflation ends when the classical back-reactions from g 2 (fp"ip 2 and A^ 4 
interactions induce large masses for 4> and ip such that they roll rapidly to the global minimum. 
Therefore, in the analysis below, we work in the limit where < 1 and consider the end of 
inflation when 6 = 6f ~ 1. 

To calculate the evolution of curvature perturbation from Eq. (|44p we need to estimate 9' and 
<5sk- The derivative of 6 in field space is 

* = (45) 

Since r ~ 1/N e <C 1, the first term is much larger than the second term by a factor of 1/r and 
6' ~ ip"/<f)'. Furthermore, as mentioned before, 5s^ = $1(1^)5^ for 9, 0' -C 1 where f2(/c) is given 
by Eq. ([33]) . Combining the above expressions for 9' and 5s, the final curvature perturbation is 
given by 

K(n f ) ~ TZo - 2 ^gVv • (46) 

There are some comments in order before we move forward. First, one can check that the integrand 
above scales like e 2/3n//3r which is fast growing so one can safely ignore the contribution from the 
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lower limit of the integral. Second, almost all /c-dependence in the expression above comes from 
Q(k). This means that the main contribution to the spectral index is induced from the entropy 
perturbations. In other words, the spectral index of the curvature perturbations and the entropy 
perturbations are more or less the same. 

Now we proceed to approximately evaluate the integral in Eq. (|46p . For this purpose note 
that (J)' scales like exp(rn) which is nearly constant so it can be taken out of the integral. Also we 
can use the background ip equation, Eq. (|15p . to replace tp" in favors of ip and ip'. After these 
simplifications, one obtains 



ml , m 2fi(fc) 



n f f n f 

+ j3 (e~ 2rn f - l) / ij) 2 &n 



(47) 



in which the function j3 (e _2rn — l) is taken out of the integral by the same reasoning as for (fJ . The 
first integral is a total derivative which can be calculated easily. To calculate the second integral, 
note that from Eqs. [T8l and TA8I one has 

-l 

n \ 

2ru ~ 



i>~ i>' (-3/2 + ru- ^ ^ ) . (48) 



This can be used to transform the second integral above into an approximate total derivative 
containing ipip' . With these simplifications employed, one obtains 

P (e~ 2rn ~ 1) f V> 2 dn * - ^3 + f (1 - e" 2 ™)^) (49) 

Plugging this into Eq. (|4T|) yields our analytic formula for the curvature perturbation 

K{n) =K, + ^{1- e~ 2 ™) n(fc)^ , (50) 

where Q(k) is given in Eq. (|33|) . 

In Fig. [T] we have plotted the predictions of the curvature perturbations from Eq. (150p and 
compared them with the full numerical results. As can be seen, they are in good agreements. 
Furthermore, the induced curvature perturbations from the entropy perturbations dominates by 
about two orders of magnitudes over the initial adiabatic perturbations. Below we find an analytic 
expression for this enhancement factor. 

Now the important question is "what is the final amplitude of the curvature perturbation?". 
As long as we are not concerned about the curvature perturbations' ^-dependence, we can find a 
simple answer for this question. For this purpose we also need to determine when inflation ends 
and the curvature perturbations saturate. So far we were working in the limit where the classical 
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N 

FIG. 1: Here we plot In \Rk\ for different modes. The blue solid curves are obtained from the full numerical 
analysis whereas the red dashed curves are obtained from our analytical formula, Eq. (|50j) . setting Of = .85. 
The induced curvature perturbations from the entropy perturbations dominate over the initial adiabatic 
curvature perturbations by about two orders of magnitudes. Note that the apparent singularity at N ~ 30 
is due to the fact that Rk vanishes at this point so In \ Rk\ diverges, otherwise it has no physically significant 
meaning. From top to bottom, the curves correspond to modes which leave the horizon at TV = 3, N = 6 
and N = 9 e-foldings respectively. The waterfall phase transition happens at N c = 7. The parameters are 
M = 7.8 x 10~ 7 to p , to = 2.5 x 10~ 7 to p and g 2 = 2A = 2 x 10~ 12 . 



back-reactions are subdominant and 6 <C 1. Once the back-reactions become important we expect 
that 6 to increase significantly. Specifically, from Eq. (|45p one observes that 

* ~ £ I (51) 

~ 0' l + tan 2 (0) ' lD ' 

so once increases significantly 9' vanishes quickly indicating that both fields approaching to their 

minima. This suggests that the time when the curvature perturbations saturate, which is nearly 

the time of end of inflation, is when 

K-'^T&H 1 - (52) 

Imposing this criteria in Eq. (|50[) and using the approximations <f>'(nf) ~ e~ rn f <f)'(n*) and ip' ~ 

77 (l — e _2rra ) ip (derived from Eq. ([48 p ) yield the following result for the amplitude of the curvature 
3 
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perturbations at the end of inflation 



(53) 

This is an interesting result. This indicates that the induced curvature perturbations from the 
entropy perturbations dominates over the initial adiabatic curvature perturbations TZq by the factor 

ip (fl f) 

— —i—e rn f 3> 1. This enhancement can be seen in Fig. [TJ Note however that this enhancing 

Wi 

factor can not be arbitrarily large. As can be seen from Eq. (|39|) the quantum back-reactions can 
become important should we start with arbitrarily small tpi. It worth mentioning that Eq. (|53p 
has no precise fc-dependence and if one is interested in ^-dependence of the curvature perturbation 
one should use the original formula Eq. (|50p with the /^-dependence dictated by O(fc). 

In Fig. [2] we have plotted the spectral index of curvature perturbations, n-ji — 1, obtained from 
our analytical formula, Eq. (|50p . compared with the full numerical analysis. As can be seen they 
are in good agreements. Both curves indicate the running from a blue spectrum to a red spectrum, 
depending on whether the mode of interest leaves the horizon before the phase transition or after 
the phase transition. This is in light of discussions below Eq. (|3ip . 



1 + 



VI. SN FORMALISM AND NON-GAUSSIANITY 



In the previous sections we have employed the perturbative approaches in details to calculate 
the curvature perturbations induced from the entropy perturbations. In this section we use the 
complementary 5N formalism to verify the previous results for the curvature perturbations. This 



also enables us to ca" 
In 5 N formalism 



cu 
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ate the non-Gaussianity parameter /tvl hi our model directly. 



5fJ] curvature perturbation on super-horizon scales can be determined by 



the variation of number of e-folds with respect to the field values at the time of horizon crossing 

C ~ n ^ 5N , (54) 

in which 

5N = N (4> + 8<f>, i> + 5ip) - N$, 4>) . (55) 

Here N, <f> and if) respectively are the background number of e-foldings and the background fields 
starting from the initial flat hyper-surface to the final constant energy density hyper-surface. One 
should note that in 5N formalism the variation should be performed with respect to the field 
values at the horizon crossing. It is also worth noting that in our case the hyper-surface of end of 
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FIG. 2: The spectral index of curvature perturbations, iir — 1, for modes which leave the horizon at 
e-folding N. The blue dotted curve is obtained from the full numerical analysis whereas the dashed green 
curve and the thin solid red curve, respectively, are obtained from our perturbative analysis, Eq. (|50p . and 
the SN formalism, Eq. (1521 . The parameters are as in Fig. [T] As explained in the text, for modes which 
leave the horizon approximately before the phase transition (N c — 7) one has a blue spectrum whereas for 
modes which leave the horizon approximately after the phase transition they become red-tilted. 



inflation is nearly the same as the surface of constant energy density. This is because the former 

aH 2 

hyper-surface is given by the relation tp 2 = — 5— where the back-reaction of the waterfall field on 



9 2 



the inflaton field terminates the slow-roll condition [35J]. From our previous analysis we know that 
the main part of energy at the end of inflation is due to ip field so a constant ip hyper-surface nearly 
coincides with a constant energy density hyper-surface. 

To employ the SN formalism we have to use the background classical trajectory. Using Eq. 
( (I18p ) one can obtain an implicit function of n in terms of fields at the time of horizon crossing 

VKn) =± V* e~^ n ~ n *^ 2 <- , (56) 

Y v (&e 



in which is the value of ip at horizon crossing, ra* . The variation of the above equation results 
in: 

= ^n_ 1 1__ 
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where 



v.; I ^ e ~ Tn 



f(n) = | + yfjie-"—}— L- . 

Y v ( &e~™ 



Here the prime denotes derivative with respect to the argument of the Bessel function. At leading 
order the curvature perturbation at the end of inflation is given by 

n ~ p^5^ , (59) 

where Sip* is the value of quantum fluctuations at horizon crossing: 5if)* = ( —== J . The power 
spectrum of curvature perturbation can be obtained by the above equations Viz = ■^j(5^ if 5'tp*) 
which results in 

(60) 

where n e = N e — N c . In our case n e is about 55 or so. This expression for the curvature perturbation 
should be compared with Eq. (|50p for the curvature perturbations obtained from the perturbative 
approach. We have checked that they are in good agreement with themselves and with the full 
numerical solution. From the above expression one can observe that the leading /c-dependence of 
curvature perturbation comes from 1/V>* factor, so the spectral index of curvature perturbation 
can be obtained as: 

dlnV* dlnV>* , . 

n * " 1 * - 2 ^inF = - 2 ^T • (61) 

Using Eq. (I57jl the final result for spectral index is the following 

y' ( v2 e -™* 

n n * 1 - 2/K) = I + ^e~ rn *^)^- f ^ (02) 

2 nf^e- 



This expression for the spectral index of curvature perturbations should be compared with Eq. 
(|3ip . the spectral index of the entropy perturbations obtained from the perturbative approach. 
As demonstrated in previous section, the main source of curvature perturbations comes form the 
entropy perturbations so the spectral index of curvature perturbation has the same form as that 
of the entropy perturbations. Furthermore, from Eq. (|62p one observes that for modes which leave 
the horizon before the phase transition, that is n* < 0, the spectral index is blue-tilted. However, 
for modes which leave the horizon (slightly) after the phase transition the spectral index is red- 
tilted. In Fig. [2] we have plotted the spectral index obtained from Eq. ([62]) and compared it with 
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10 20 30 40 50 



N 

FIG. 3: Here we present ln|/jvx| for different modes obtained from the full numerical analysis. The solid 
blue curve, the dashed red curve and the dashed-dotted green curve, respectively, correspond to modes 
which leave the horizon at N = 3, N = 6 and N = 9 c- foldings. As can be seen, the final value of Jnl 
saturates at $nl — 0.2 which is very well approximated by our analytical formula Eq. (|64[) . The bump at 
N ~ 30 is because Rk = at this point and since Jnl ~ RZ 1 so it peaks at this point. Otherwise, it has 
no physically significant meaning. Similarly, the apparent singularity at iV~40 is because /jvt changes sign 
at the final stage of inflation and ln|/jvx| diverges at this point. As before, it has no significant meaning 
because we measure fwL and Rk at the end of inflation. The parameters are the same as in Fig. [T] 



the spectral index obtained from the full numerical analysis and from the perturbative analysis, 
Eq. (|50p. The running from a blue spectrum to a red spectrum is common in all three curves. 

Finally we can easily obtain the non-Gaussianity parameter in our model by the following 
formula: 



6 

TJNL = — 
5 n 



fNL = ^f^ ~f(n e ) . (63) 



Using the approximate behavior of Bessel functions given in the appendix, Eqs. IA7I and (|A8j) . to 
simplify f(n e ) we obtain 

t 5 ( 3\ 5/3 

TNL — ~ I TV C± . 64 

6 V 2 J 18 v ' 
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This is a very interesting formula indicating that the level of non-Gaussianity in our model is 
controlled by the parameter (3 measuring the sharpness of the phase transition during inflation. 
For our numerical example with (3 = 0.7, we obtain /jvx — 0.2. In Fig. [3] we have plotted /jvx 
obtained from our full numerical analysis for different modes. As we see, /jvl measured at the end 
of inflation is very well approximated by our analytical estimation Eq. (|64|) and saturates at the 
expected value /nl = 0.2. Also as can be seen from the plot, /nl has no significant fc-dependence. 

We see that generically /jvx < 1 in our model because of the assumption that the phase 
transition is mild during inflation. However, this level of non-Gaussianity is much bigger than non- 
Gaussianity predicted in simple inflationary models which is at the order of the slow-roll parameters. 
Increasing j3 by one or two orders of magnitude one can obtain significant non-Gaussianities 
However, with /3 > 1 we approach the standard hybrid inflation model where the waterfall phase 
transition is very sharp and the quantum back-reactions dominate exponentially over the classical 
back-reactions and one should take their effects into account. This conclusion is also consistent 
with the bound obtained in Eq. ()42|) in order to neglect the quantum back-reactions compared to 
the classical back-reactions. 



VII. CONCLUSION AND DISCUSSIONS 



In this work we considered a variant of hybrid inflation where the waterfall phase transition is 
mild such that a long period of inflation can be obtained after the phase transition. We found the 
model parameters where the quantum back-reactions are sub-dominant compared to the classical 
back-reactions. As can be seen from Eq. (|42|) the strength of quantum back-reactions is exponen- 
tially sensitive to the parameter /3. For a fixed values of inflaton mass and the coupling g, one has 
to start with small enough f3 such that the quantum-back-reactions are negligible. 

We have calculated the curvature perturbations and the spectral index perturbatively as well as 
using the complementary 5N formalism. We have shown that the curvature perturbations induced 
from the entropy perturbations dominate over the adiabatic perturbations by the enhancing factor 
&e rrif 1. Consequently, the spectral index of curvature perturbations has the same form as the 
spectral index of the entropy perturbations. We have shown that n-ji runs from a blue spectrum 
to a red spectrum depending on the time when the mode of interest leaves the horizon. For the 
modes which leave the horizon approximately before the phase transition n-n > 1 whereas for the 
modes which leave the horizon approximately after the phase transition n-ji < 1. This may have 
interesting consequences when compared to the WMAP data [1] to fit the CMB observations. We 
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would like to come back to this question in future. 

Using the <5iV formalism we have calculated the level of non-Gaussianity in our model. We 
found the interesting result that /jvl — ff ^ 1- This indicates that the parameter [3, which is a 
measure of the sharpness of the phase transition, controls not only the strength of the quantum 
back-reactions but also the level of non-Gaussianity. Although this level of non-Gaussianity may 
not be detectable observationally in near future, but it is one order of magnitude larger than the 
level of non-Gaussianity predicted in simple models of inflation. 

As studied in [30j], in our model with (3 < 1, primordial black holes can be formed copiously 
which may be dangerous cosmologically. The analysis of black hole formation and their cosmological 
consequences are beyond the scope of this work. We would like to come back to this question in 
future. 
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Appendix A: Approximate behaviors of Bessel functions 

In this appendix we present some asymptotic behaviors of Bessel functions for the range of 
parameters relevant to our analysis. These approximate relations are very useful in evaluating C2 
in Eq. (j27|) which is also used to calculate the curvature perturbations power spectrum and the 
spectral index. 

The Bessel functions J u (x) and Y u (x) satisfy the following equation of motion 



We are in the limit where y > i > 1. In this limit the above equation simplifies approximately to 




(Al) 



/"(*) - -f(x) ~ 



(A2) 
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which has a simple solution of 

/(*) ~ cx ±(v / 4^TT+l)/2 „ cx ±u (A3) 

Here c is a constant of integration and throughout this appendix the upper and lower signs belong 
to the Bessel function of the first and second kind, respectively. To improve our approximation, 
note that the main error in the analysis above is in ignoring the factor 1 in comparison with v 2 jx 2 
in (lAip . Therefore, to next leading order, the differential equation is 

/''(*) +(l-^W) = 0. (A4) 

This equation modifies the solution ()A3|) by allowing c to be x-dependent, whereas previously it 
was a constant. Since this is the second order approximation we can assume that c(x) is a slowly 
varying function of x and ignore its second derivative. Using this assumption and substituting 
(|A3P in to (|A4j) results in the following equation for c, 



iyV 2 + 1 + 1) c'{x) + xc(x) ~ 0, (A5) 



x 2 \ „ ( x 2 



which has a simple solution 

c(x) ~ cexp I ; ) ~ cexp (——-). (A6) 

V ' V ±2\Z4i/ 2 + l + 2/ V K ' 

Finally, we should determine the constant c. Since the overall behavior of our approximate solution 

looks like the standard small argument limit of Bessel function (i/~i<1), this can be used to 

fix c so our approximate solution conforms the small argument expansion of the Bessel function. 

This leads to our desired formulae 

Mx) = f(^M)(f)" exp ("^ 

(!)"-(£)' (A7) 

One can check numerically that the above equations are good approximation of Bessel functions. 

Now we obtain an approximate expression for c\ and C2 in Eq. ( (|27p ) which we need for the 
calculation of spectral index and quantum fluctuations. From (|A7p one has, 



Yl{x) ~ Y v {x) + !;) , (A8) 

where the prime denotes the derivative with respect to the argument x. Using these expressions 
and also the explicit form of Hankel functions one obtains the following approximations for c\ and 
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(>2 



C1 



(>2 



2 ire- Na Y e 
k r 8r 1 z 



v 3/2 



■2u 2 -x 2 , 
-r(i + z) I + 3z + 3i 



2 vre" 7 ^ e iz 



ru + z) 



2i/ 2 - x 2 
v 



+ 3z + 3i 



(A9) 



Here x and z are the arguments of Bessel and Hankel functions in Eq. ( (|27p ). respectively, denned 
by 

_ _ VA)(r + l)„-rn m 



— e~ nm ~ v / 2^or < 1 . 



(A10) 
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